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Strain-Dependence of Surface Diffusion: Ag on Ag(lll) and Pt (111) 
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Using density-functional theory with the local-density approximation and the generalized gradient approx- 
imation we compute the energy barriers for surface diffusion for Ag on Pt (111), Ag on one monolayer of Ag 
on Pt (111), and Ag on Ag (111). The diffusion barrier for Ag on Ag (111) is found to increase linearly with 
increasing lattice constant. We also discuss the reconstruction that has been found experimentally when 
two Ag layers are deposited on Pt(lll). Our calculations explain why this strain driven reconstruction 
occurs only after two Ag layers have been deposited. 
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q ■ The diffusion of a single adatom on a surface is a funda- 

mental process determining the surface morphology de- 
veloping in epitaxial growth. For Ag on Ag(lll) it has 
been demonstrated!!! that the substrate temperature can 
be lowered periodically for a short period of time so that 
diffusion is suppressed, and as a result the density of is- 
lands is increased at the beginning of growth of a new 
layer. With this trick the growth mode can be changed 
from three-dimensional to two-dimensional (for a discus- 
sion of this apd other approaches affecting the growth 
mode see RefBI). Growth of one material on a different 
material is of particular interest for a number of tech- 
nological applications. In such a heteroepitaxial system 
the material to be deposited is under the influence of epi- 
taxial strain, yet very little is known about the influence 
of strain on the surface diffusion constant. Additionally, 
strain is not only due to lattice mismatch but it is also 
present on surfaces of homoepitaxial systems as a result 
of the modified bonding configuration. 

Experiments show that the diffusion constant follows 
an Arrhenius type behavior D ~ exp(— Eb/ksT), where 
Ef, is the barrier for surface diffusion and T is the sub- 
strate temperature. In a recent scanning, tunneling mi- 
croscopy (STM) experiment Brune et alEl measured the 
island density N of Ag on Ag(lll), Ag on Pt (111), and 
Ag on one monolayer (ML) of Ag on Pt Q 11) for different 
temperatures. With the scaling relatione! N ~ (D/F)~ x , 
where F is the deposition flux, the diffusion barrier Eb 
can be n determined as long as the temperature is low 
enoughO so that the scaling exponent \ = 1/3- Brune 
et al. find that the diffusion barriers are E, s ~ = 
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157± 10 meV for Ag on Pt (111), E l 



Ag-Ae 



97±10meV 



,Ag-Ag/Pt 



for Agon Ag (111), and Ep ^ ,rh = 60±10 meV for Ag 
on one ML Ag on Pt(lll). Thus, the diffusion barrier 
for Ag on top of a pseudomorphic layer of Ag on Pt (111) 
is substantially lawer than that for Ag on Ag(lll). It 
has been arguedtl that a metallic monolayer which is 
supported on a dissimilar substrate is electronically per- 
turbed by the substrate, and that its chemical properties 
are altered. The question arises whether the low diffu- 
sion barrier for Ag on one ML Ag on Pt (111) is a result 



of the compressive strain of 4.2% or an electronic effect 
due to the Pt underneath the Ag layer. 

There are only a few theoretical studies that inves- 
tigated the effect of lattice mismatch on the diffusion 
barrier. In an extensive molecular dynamics study for 
Si on Si (001) that employed a Stillinger- Weber potential 
Roland and GilmerQ found that the barrier for diffusion 
along the fast channel parallel to the dimer rows is low- 
ered by approximately 10% for both, 3% compressive 
or 2 % tensile strain, while diffusion along the same di- 
rection atop the dimer rows is increased by about 10 %, 
so that a general trend cannot be seen. Using Lennard- 
Jones potentials Schroeder and Wolf! found that the dif- 
fusion barrier increases linearly as the lattice constant 
increases. For a metallic system we are only awaxe of re- 
sults for Ag on Ag (111) where the authors of Refill find in 
an effective medium theory (EMT) calculation that the 
diffusion barrier increases with increasing tensile strain 
of the surface and decreases with increasing compressive 
strain. Comparison with experiment [for Ag on Pt (111), 
Ag on Ag (111), and Ag on one ML Ag on Pt (111)] con- 
firms this trend but also shows that the EMT results are 
off with an error between 20% and 100%. For values of 
misfit larger than 3 % the EMT diffusion barrier starts 
to decrease, which is in qualitative disagreement to our 
results described below. 

In this work we present first-principles calculations of 
the dependence of the diffusion barrier on the lattice con- 
stant for Ag on Ag(lll). This allows us to isolate the 
effect of strain from electronic effects caused by the dif- 
ference of materials. The dependence of the diffusion 
barrier on lattice mismatch is found to be essentially lin- 
ear. The diffusion barriers for the system Ag on Pt (111) 

are calculated as well. We find E b s ~ = 150 meV, 
^Ag-Ag = gl meV; and ^Ag-Ag/Pt = 63 meV ^ h . ch . g 

in good agreement with the experimental datau within 
the error margins. 

We employ density- functional theory (DFTl. together 
with the local-density approximation (LDA)t3 for the 
exchange-correlation (XC) functional. The energy bar- 
riers are found to be only weakly affected when the 
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FIG. 1. A schematic representation of the adsorption sites 
(fee and hep site) and saddle point (bridge site) on the (111) 
surface. 

generalized-gradient approximation (GGA)E9 is used in- 
stead of the LDA. In this study the GGA results for 
the barriers are typically higher but not more than 
5 % to 10 % compared to the LDA results. Norm- 
conserving, fully-separable pseudopotentials have been 
employed that were generated according to a scheme pro- 
posed by Troullier and Martins .til Thej-computer code 
used is described by Bockstcdte et alil3. The calcu- 
lated bulk lattice constants are apt = 3.92 A for Pt and 
a Ag = 4.05 A for Ag with the LDA and a Pt = 4.01 A 
and OA g = 4.19 A with the GGA. The LDA results axe 
smaller than the GGA results which is the usual trend.E-3 
For Ag the experimental value is slightly larger (smaller) 
than the LDA (GGA) result while for Pt the LDA re- 
sult is identical to the experimental lattice constant. In 
the above theoretical values the influence of zero-point 
vibrations is not included; it would increase the lattice 
constant by less than 0.2 %. 

To simulate the surface we use the supercell approach. 
In the z-direction slabs are separated by a vacuum region 
with a thickness that is equivalent to 6 layers, and it has 
been tested carefully that this vacuum region is thick 
enough. The adatom is placed on only one side of the 
slab. Because of the simple geometry of the (111) surface 
(cf. Fig. 1) 

only the fee and hep sites need to be considered as 
adsorption sites with the bridge site as the saddle point. 
It turns out that the fee site is always slightly favored 
over the hep site so that all barriers discussed are the 
differences between the total energies of the system with 
the adatom in the fee site and the bridge site. 

We carried out careful tests varying the slab thickness 
and surface cell size to ensure that interactions between 
neighboring adatoms are negligible. The tests for Ag on 
Ag(lll) are summarized in Table 1. In calculation 1, 2, 
and 3 we varied the slab thickness Mayer and conclude 
that M a y Cr = 4 is sufficient. The electronic wave func- 
tions are expanded in plane-waves that are truncated at 
a cut-off energy E cu t. From calculation 2, 4, and 5 we 
find that E cut = 40 Ryd is sufficient if the desired accu- 
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FIG. 2. The diffusion barrier as a function of the relative 
lattice constant a/arj for Ag on Ag (111). The lattice constant 
is normalized to the computed bulk lattice constants of Ag 
a = 4.05 A(DFT) and a = 4.075 A(EMT). The EMT results 
are taken from Ref. 3. The solid lines are a guide to the eye. 



racy is ±10 %, but for a accuracy of ±2 % a larger cut-off 
-E^cut = 50 Ryd is necessary. For all the results reported 
below we choose a slab with a (2 x 2) cell as it can be 
justified from calculation 2 and 7. We always relaxed the 
positions of the adatom and the atoms of the top layer. 
A test revealed that the results remained unaltered upon 
relaxation of the second layer. For the k-summation we 
used 10 k points in the irreducible part of the surface 
Brillouin zone of the (2 x 2) cell that were generated ac- 
cording to Ref. [l|. The k points included the T-point 
but calculations with different k points that did not in- 
clude the T-point (calculation 6) yielded identical results. 
Thus, the error in the Brillouin zone integration is negli- 
gible. 

Results for the diffusion barrier for Ag on Ag(lll) as 
a function of lattice mismatch are shown in Fig. 2. For 



the unstrained system we obtain E. 



Ag-Ag 



81 meV in 



good agreement with the STM results of Ref. |[ In the 
range of the lattice constants studied the diffusion barrie* 
Eb(a/ ao) varies linearly within our numerical accuracy^ 
with a slope of ~ 0.7 eV. This is in contrast to the EMT 



TABLE I. Convergence Tests for 
fully self-consistent GGA and GGA- 
described in the text. The cell size 
in calculation 7 which was obtained 
the number of k points. 



Agon Ag(lll). GGA is 
ap is a posterior GGA as 
used was (2x2) except 
with a (3 x 3) cell A^ is 
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FIG. 3. The total energy of an adatom in the fee site and 
bridge site as a function of the relative lattice constant a/ao 
for Ag on Ag (111). The lattice constant is normalized to the 
computed bulk lattice constants of Ag ao = 4.05 A (DFT) 
and a = 4.075 A (EMT). The solid lines are a guide to the 
eye. 



results of Ref. |3j that are also shown in Fig. 2 where 
the dependence is sub-linear and the diffusion barrier de- 
creases for misfits larger than ~ 3 %. 

To further explain the strain dependence of the diffu- 
sion barrier we show in Fig. 3 the strain dependencies 
of the underlying contributions, i.e., the total energies of 
an adatom at the fee hollow and the bridge sites. These 
are computed by subtracting the total energies of a clean 
surface and a free silver atom from the total energy of 
a surface with an adsorbed atom. The total energies, 
of Fig. 3 show an approximately linear dependence!!^ 
within a surprisingly large range of lattice constants. The 
change of the total energy E^ia/a®) with the adatom at 
the threefold coordinated fee site is stronger than that 
at the twofold coordinated bridge site, and the slopes 
are ~ 2.1 eV for the bridge site and ~ 2.8 eV for the 
fee site. The difference reflects that the corrugation of 
the potential energy surface gets more pronounced when 
the surface is expanded while under compression it be- 
comes more flat. Obviously, the value of the difference 
equals the value of 0.7 eV .obtained above from Fig. 2. 
We note that Dobbs et aZO recently pointed out that 
AEb = (cTf cc — rTb r idge)e where e is the strain and <7f cc 
(rjbridge) is the surface stress with and adatom at the fee 
(bridge) site. The slopes of the curves in Fig. 3 give these 
stresses <7f cc and abridge- 

The diffusion barrier for Ag on Pt (111) has been cal- 



culated as E, 



-pt 



150 meV with a (2 x 2) cell, a slab 



with 4 layers where the atoms of the top layer and the 
adatom are relaxed, and E cut = 40 Ryd. The same result 
has been obtained for i) a thickness of 3 layers, it) the 
top two layers relaxed, and Hi) a cell size of (3x3). We 
thus conclude that our result is converged. The value for 



Ag-Pt 

b 

Ag-Pt 



is in excellent agreement with the STM value of 

=|A57 ±10 meV. Previous DFT-LDA results by 

FeibelmantS who used the Green function theory gave a 



higher value of E t 



Ag-Pt 



200 meV. 



Our calculations predict that the diffusion barrier 
drops dramatically for Ag on Pt (111) after one complete 
layer of Ag has been deposited. For this diffusion barrier 
we obtain E, s s ' = 63 meV in good agreement with 
the STM value of E Ag ~ Ag/Pt = go ± 10 meV. For Ag on 
one ML of Ag on Pt (111) we have chosen E cut — 40 Ryd 
and 3 Pt layers. We also checked the result with only 2 
Pt layers and verified that choosing 3 Pt layers is indeed 
sufficient for this system. The diffusion barrier for Ag on 
one ML Ag on Pt (111) is almost identical to the value 
of E Ag ~ As = 60 meV we obtained for Ag on Ag(lll) 
compressed to the lattice constant of Pt. This suggests 
that the small diffusion barrier for Ag on one ML Ag on 
Pt (111) is mainly due to the effect of strain and not an 
electronic effect because of the Pt underneath. 

All the results discussed above have been obtained with 
the LDA for the XC functional. It is an ongoing debate 
whether the GGA is really an improvement over the LDA 
and we investigated the influence of the GGA for all the 
systems discussed here. For Ag on Ag (111) we find that 
the effect of the GGA on the diffusion barriers is only 
about 5 meV. This can be seen in table 1 from calcula- 
tions 4 and 8 for a self-consistent GGA according to Ref. 
pT| and from calculations 2 and 9 where the energy has 
been computed with the GGA from a self-consistent LDA 
electron density (GGA-ap) according to Ref. |l2|. Simi- 
larly the barrier increases by only 5 — 10 meV for Ag on 
Pt (111) and Ag on one ML Ag on Pt (111). The negli- 
gible effect of the GGA is in agreement with results by 
Boisvert et al. for self diffusion of Pt on Pt (111)E3 and 
results by Yu and Schcfflcr for diffusion via the hopping 
mechanism for Ag on Ag(100).EJ The (111) surface is a 
close packed surface with a very small surface corrugation 
and since LDA and GGA results are very close it is plau- 
sible to assume that for such a system LDA and GGA are 
both good approximations for the exact XC functional. 
This situation will presumably be different for diffusion 
events along and especially -across steps as it has been 
found by the authors of Ref.E3 for Ag on Ag (100). 

Brune et alE also measured the island densities of Ag 
on two ML of Ag on Pt (111) and found that the island 
density is much larger than it is for Ag on just one ML 
of Ag on Pt(lll). The reason for this increased island 
density is not a higher barrier for surface diffusion. The 
second layer of Ag on Pt (111) reconstructs in a trigonal 
network where-domains with atoms in the fee and hep 
site alternate.EIJ This reconstruction occurs either during 
growth with high enough adatom mobility or upon an- 
nealing and it can be concluded that this trigonal network 
is the equilibrium structure. The periodicity of these do- 
mains is approximately two domain boundaries for every 
24 atoms. This can be understood with purely geomet- 
rical arguments because the lattice mismatch is ~ 4.2% 
and every domain boundary implies that there is half of 
an Ag atom less, so that the domain network provides an 
efficient mechanism to relief epitaxial strain. The exper- 
iments indicate that these domain walls act as repulsive 



walls so that the island density is determined by the size 
of the reconstructed unit cell and not the barrier for self 
diffusion on the flat terrace. It was not clear however 
why this domain network is formed only after 2 ML Ag 
have been deposited and not already upon completion of 
the first Ag layer. 

To answer this question we compared the total energy 
of an adatom in the fee and in the hep site. Calculations 
were carried out with a (1 x 1) and a (2 x 2) cell (i.e. 
coverages 9 = 1.0 and O = 0.25) and slab thicknesses of 
up to 5 layers. We find that the fee site is energetically 
more favorable than the hep site in all cases. The energy 
difference between those two sites for an Ag adatom on 
Pt (111) is 40 meV for a (2 x 2) cell (30 meV for a (1 x 1) 
cell), but less than 10 meV for an Ag adatom on one ML 
Ag on Pt (111). The reconstruction does not occur in the 
first layer because the Ag atoms are bound much stronger 
at the fee sites but does occur in the second layer because 
the total energies for adatoms in the fee and hep site are 
almost indistinguishable. 

In conclusion we found that in the range studied here 
the diffusion barrier E^a/ao) for Ag on Ag(lll) in- 
creases approximately linear with a slope of ~ 0.7 eV 
when the lattice constant increases. We propose that 
this result might offer a possibility to change the diffu- 
sion barrier and thus the island density during growth 
by artificially straining the substrate. If this can be done 
periodically in a controlled manner it might be an alter- 
native approach to vary the diffusivity in situ. Our LDA 
results show that careful DFT calculations with the LDA 
reproduce experimental values for the surface diffusion 
barrier for the system Ag on Pt (111) with very high ac- 
curacy. Work to calculate the barriers for diffusion along 
and across steps and across domain boundaries and pre- 
factors is in progress. 
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